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ABSTRACT 



We present a scheme to determine the spectral state of the canonical black hole Cyg X- 1 using data from previous and current X-ray 
all sky monitors (RXTE-ASM, Swift-BAT, MAXI, and Fermi-GBM). State determinations of the hard/intermediate and soft state 
agree to better than 10% between different monitors, facilitating the determination of the state and its context for any observation of 
I the source, potentially over the lifetimes of different individual monitors. A separation of the hard and the intermediate state, which 

strongly differ in their spectral shape and short-term timing behavior, is only possible when monitor data in the soft X-rays (<5 keV) 

qq are available. A statistical analysis of the states confirms the different activity patterns of the source (e.g., months to years long 

hard state periods or phases during which numerous transitions occurs). It also shows the hard and soft states to be stable, with the 
probability of Cyg X-1 remaining in a given state for at least one week to be larger than 85% for the hard state and larger than 75%, 
for the soft state. Intermediate states are short lived, with a probability of 50% that the source leaves the intermediate state within 
• three days. A reliable detection of these potentially short-lived events is only possible with monitor data with a time resolution of 

CO better than 1 d. 

O 

Keywords, stars: individual: Cyg X-1- X-rays: binaries - binaries: close 

> 

• h 1. Introduction In a hardness intensity diagram (HID) transient BHBs move 

/\ on a clear trajectory, the so-called q-track dFender et al. 120041. 

^ Accreting galactic black hole binaries (BHBs) show two main Different states correS pond to different parts of the track. Of 
Cd spectral states: a soft state with a thermal X-ray spectrum dom- ial interest ifj the so . called jet . line that M coincides 

mated by an accretion disk and a hard state with a power law wkh the transition from hard interm ediate to soft intermediate 
spectrum with a photon index r ~ 1.7. The intermediate or torn- sMe and is associated with radio ejection events Render et al. 



sitional state can be subdivid ed into a hard intermediate state A q . track .i ike behavior on HIDs or on their generaliza 
and a soft intermediate state (Bellom 201 01. Transient B HBs \^ the disk-fraction luminosity diagrams, has been found in 
also show a quiescent state ( |McClintock & Rirrul lard 2006). All different accre ting sources such as neutron star X-ray binaries 
states show distinct spectral and timing properties in the X-rays. ( iMaitra & Bailyn||2004| Aql X-1), dwarf novae flKordingl 
Radio emission is detected in the hard and hard inter mediate | et al.|2008 SS Cyg), or AGN (iKording et al.|2006|. It is there- 
states and i s strongly suppressed in the soft states (e.g., |Fender| ' fore kely to reflect basic accretion/ejection physics inherent to 
|et al.|2009|>. a w j de ran g e D f accreting objects. 

Send offprint requests to: V. Grinberg, Cygnus X-1 is a key source for understanding accretion and 
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relatively nearby with a radio parallax distance of 1.86^q'Jj kpc 
( Reid et al.|201 1 consistent with X-ray dust scattering halo es- 
timates of Xiang et al.|201 1[ ), bright (above ~100mCrab in the 
1.5-12 keV band in the hard state), and often undergoes (failed) 
state transitions (e.g., Pottschmidt et al. 2003 1, which are thought 
to be connected to changes in its radio jet (Fender et al. 2004, 
[20061 |Wilms et al.1[2007l i. The bl ack hole is in a 5.6 d o rbit 
around its donor star, HDE 226868 (Brocks opp et al.|1999b| and 
references therein). The orbital modulation is also detected in 
radio (e.g., |Pooley et al.|T 999) and, due to modulation of the soft 
X-ray flux by absorption in the donor's stellar wind, in the X- 
rays (Balucihs ka-Church et al.|2000||Poutanen et al.|2008|>. Ad- 
ditionally, CygX-1 shows a superorbital per iod of 150 d (|Brock-| 
sopp et al.||1999a| |Benlloch et al.||2004} |Poutanen et al.||2008) 



although that superorbital variabi 



ity seems to be unstable and 



has recently been reported to have doubled to 300 d (Zdziarski| 



et al. 201 1 1. In the hard state, radio jets have been observed (e.g., 
Stirling et al.|200l| ). The emisssion above ~400keV CygX-1 is 



strongly polarized ( |Laurent et al.|201 l]|Jourdain et al.|2 012). 



As a persistent source CygX-1 does not cover the full q- 
track: its bolometric luminosity changes only by a factor of ~3~ 
4 between the states ( [Wilms et al.|2006] and references therein) 
and its spectrum is never fully disk-dominated. The frequent 
state transitions (sometimes very fast - within hours, see, Bock 
|et al.|20lT) mean that the source often crosses the jet line. Since 
these state changes are thought to be associated with significant 
changes in the accretion flow geometry and energetics, a knowl- 
edge of the source state is crucial for the interpretation of all 
(multiwavelength) observations of CygX-1 and its donor star. A 
typical example is the study of the stellar wind of HDE 226868, 
which during soft states is strongly photoionized by the radiation 
from the vicinity of the black hole ( |Gies et al.|2008) >. 



In the past decade, state information was readily available 
using the All Sky Monitor on the Rossi X-ray Timing Ex- 
plorer (RXTE-ASM) and regular pointed monitoring observa- 
tions. Here, various state definitions exist, which use, e.g., mea- 
sured count rates and/or colors (e.g., |Remillard|2005| Gies et al. 
2008), or sophisticated mapping between these measurements 
and spectral parameters (e.g., Ibragi mov et al.|[2007| |Zdz iarski 
et al.|20rT) . The former prescription is easy to use, but is very 



instrument specific and cannot be translated easily to other X- 
ray all sky monitors. The latter approach requires a sophisti- 
cated knowledge of the instrumentation of all sky monitors as 
well as of the detailed spectral modeling. Furthermore, the pre- 
viously used state definitions are all slightly inconsistent with 
each other. In this work we introduce a novel approach to clas- 
sify states of CygX-1 using the all sky monitors RXTE-ASM, 
MAXI, Swift-BAT and Femi-GBM based on 16 years of pointed 
RXTE observations. Our aim is to find an easy-to-use prescrip- 
tion for the determination of states that is as consistent as pos- 
sible between these instruments in order to facilitate long-term 
studies that exceed in duration the lifetime of individual moni- 
tors. We start with a description of our data reduction approach 
in Sect. [2] Section [3] comprises the actual state mapping from 
pointed RXTE observations to RXTE-ASM, Swift-BAT, MAXI 
and Fermi-GBM, including a discussion of the precision of state 
determinations attainable with these instruments. We summa- 
rize and discuss our results in the light of the statistics of the 
state behavior of CygX-1 in Sects. |4]and|5] 



2. Observations and Data Analysis 

2.1. ASM data 

The RXTE-ASM instrument consisted of three Scanning 
Shadow Cameras (SSCs), in which a position sensitive propor- 
tional counter was illuminated through a slit mask ( |Levine et al.| 
1996). A typical source was observed at randomly distributed 
times 5 to 10 times a day in three energy bands roughly corre- 
sponding to 1.5-3.0 keV (band A), 3.0- 5.0keV (band B), and 
5.0-12keV (band C) ( |Levine et al.fT996T l. 

We consider a ll 97556 RXTE All Sky Monitor (ASM, 
|Levine et al.|1996] l measurements of Cyg X-l performed during 
the lifetime of RXTE. In the following all data analysis was per- 
formed with ISIS 1.6.2 ( |Houck & Denicola|2000[ |Houck|2002] 
INoble & Nowak|2008D . 

We first filter for measurements where the background was 
clearly oversubtracted (count rates r^, rg, or yq < 0). The qual- 
ity of the ASM data stalled deteriorating after about MJD 55200 
(early 2010 January): valid pointings are fewer and the overall 
variance of the measurements becomes larger (Fig. [TJ. Since 
the larger variance could be source intrinsic, we analyze the 
ASM light curve of a source known to be roughly constant at the 
level of precision required here, the Crab (Wilson-Hodg e~et~aL| 
2011): after MJD 55200 the Crab light curve shows prolonged 



gaps (Fig. |2j. Where data exist, the average values of the ASM 
count rate (calculated on a 10 d timescale), which were stable at 
~75 cps before, decrease by up to 10%. The variance of the data 
strongly increases by a factor of ~10. The ASM hardness (be- 
tween the C and A bands) increases from stable values around 
~0.95 to up to ~1.4 and shows higher variability. The timescale 
of this change in behavior as well as the similar behavior of the 
light curves of other sources during the same time period imply 
that we do not see the long term variability of the Crab nebula as 
has been observed by Wilson-Hodge et aT7| ( |201 1) > here, but truly 
instrumental decline during the last years of the ASM lifetimeF] 
The ASM data measured post 2010 January can still be used 
to assess trends with large changes in count rate; for the fol- 
lowing analysis, however, which relies on absolute values of 
both count rate and hardness, we ignore all ASM data after 
MJD 55200. Overall, we use 94068 individual ASM measure- 
ments corresponding to 80082 individual times (a source can be 
in the field of view of two SSCs simultaneously) spanning over 
5000 days from MJD 50087 (1996 January 5) to MJD 55200 
(2010 January 4). 

2.2. Pointed RXTE observations 

We consider all pointed RXTE observations of Cyg X-l made 
during the RXTE lifetime (MJD 50071 to MJD 55931) Data 
from the Proportional Counter Array (PCA, Jahoda et al.{2 006) 
and the High-Energy X-Ray Timing Experiment (HEXTE, ]Roth-] 
|schild et al.l[T998l ) onboard RXTE were reduced with HEA- 
SOFT 6.11 following the procedure introduced in the previous 
papers of this series (Pottsch midt et al.| 2003 ; Wil ms et al.|2006) . 

For spectral analysis of the PCA data, we use data from 
the top xenon layer of Proportional Counter Unit (PCU) 2 only. 
PCU 2 is the best calibrated PCU and had been running during 
all RXTE observations of Cyg X-l. We make use of the im- 
proved PCA background models and therefore only discard data 



within 10 minutes of the South Atlantic Anomaly (SAA; Fiirst 



1 Vrtilek & Boroson 


l2013i come to similar conclusions regarding 


ASM data from 2010 


on; they also report on gain changes in the last 


two years of ASM lifetime as a possible cause. 
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Fig. 1. Pointed RXTE observations and all light curves (RXTE- ASM, MAXI, Swift-BAT and GBM) of CygX-1 used in this analysis. ASM, 
MAXI and BAT data are shown in the highest available resolution, GBM are binned daily. ASM hardness is calculated by dividing count rates 
in band C (5.0-12keV) by count rates in band A (1.5-3.0 keV). Vertical dashed lines and horizontal arrows represent periods of different source 
activity patterns. Blue, green, and red colors represent states of individual measurements classified using the respective classification for the 
different instruments as introduced in Sect. [3] blue represents the hard state, green the intermediate state and red the soft state. ASM data after 
MJD 55200 (shown in grey) are affected by intrumental decline. Hard and intermediate states cannot be separated in BAT and GBM; BAT and 
GBM data corresponding to these periods of hard or intermediate states are therefore shown shown in black. 
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Fig. 2. Crab ASM light curve (left panel, 1.5-12keV, gray points) and hardness (right panel, 5.0-12keV/1.5-3.0keV, gray points) with the 
respective 10 day average values (black histograms) and variances (blue histograms). 



Article number, page 3 of 12 



|et al.|2009] > passages, as opposed to the 30 minutes we used pre- 
viously. 

No HEXTE spectra are available during some observations 
either due to the non-standard observation mode employed, 
mainly in early observations which are not part of our bi-weekly 
campaign, or due to the failure of the rocking mechanism of the 
HEXTE clusters in the last years of the RXTE lifetime. 

We extract all PC A spectra in the standard2£ mode for 
each RXTE orbit and obtain 2741 individual spectra. The im- 
proved PCA background and response allow us to consider a 
wider energy range than previously, namely ~2.8 keV to 50keV. 
The wider overlap between the PCA and HEXTE spectra, which 
are considered between 18keV and 250 keV, leads to a better 
constraint on the multiplicative constant which accounts for the 
differences in the flux calibration of the two instruments. Fol- 
lowing B ock et al"1 ( 201 1) and Hanke (201 lit we add a systematic 
error of 1% to the fourth PCA bin (2.8-3.2 keV) and of 0.5% to 
the fifth PCA bin (3.2-3.6 keV) in the data from PCAs calibra- 
tion epoch 5 (2000 May 13-2012 January 5f\ which constitute 
the majority of our dataset. For epochs 1 to 4 we aim for the 
closest possible match in energy for the systematic errors. 

Where possible due to the available PCA modes, light curves 
with a time resolution of 2~ 9 s [~2 ms) were extracted for timing 
analysis, choosing the same energy bands as |Bock et al.| ( |20lT| ): 
a low energy band corresponding to energies 4.5-5.7 keV (chan- 
nels 11-13) and a high energy band corresponding to energies 
9.4-14.8 keV (channels 23-35; energy conversion for epoch 5). 
The calculation of the rms and X-ray time lags follows [Nowak| 
|et al.| ( fl999"l ). We use segments of 4096 bins, i.e., we calculate the 
root mean square variability (rms) between 0.125 and 256 Hz. 

For two simultaneous, correlated light curves, such as the 
high and low energy lightcurves of CygX-1 used here, one 
can calculate a Fourier-frequency dependent time lag between 
the two from the Fourier phase lag at any given frequency 
( |Nowak et al.||1999] and references therein). In our calcula- 
tions a positive time lag means that the lightcurve in the hard 
energy band is lagging the soft. The timelag depends strongly 
on Fourier-frequency and shows a complex variation with state 
(see [Pottschmi dt et al.|2000 for examples and note that we use 
the same frequency binning in the relevant frequency-range). To 
obtain a single value which serves as a good signature for the 
overall level of the lags, we average it over the frequency range 
of 3.2-10 Hz (see |Pottschmidt et al.|2000l [20031 . 



2.3. Swift-BAT, MAXI, and Fermi-GBM data 

RXTE was switched off on 2012 January 5 (MJD 55931). To 
continue the monitoring of the long term behavior of Cyg X- 1 
we therefore need to use other instruments. The all sky monitor 
available in the soft X-ray band at the time of writing is MAXI, 
the hard X-rays above lOkeV are covered by Swift-BAT and the 
Fermi-GBM. 

MAXI is an all sky monitor onboard the Japanese module 
of the International Space Station (Matsuoka et al.|2009| >. Light 
Curves from the Gas Slit Camera detector (GSC) are available 
in three energy bands (2-4 keV, 4-10keV and 10-20 keV) on a 

2 Over its lifetime the PCA saw four different gain calibration epochs, 
followed by a long fifth epoch defined by the loss of the propane layer in 
PCU0 in 2000 May. See http://heasarc.gsfc.nasa.gov/docs/ 
xte/e-c_table .html for details. These epochs are not to be confused 
with the five activity periods of CygX-1 that we define in Fig. [T] and 
Sect.[3H 



dedicated websit^] MAXI light curves show prolonged gaps of 
several days due to observational constraints. 



Swift-BAT is sensitive in the 15-150 keV regime (Barfhelmy 
|et al.|[2005] l. Regularly updated, satellite-orbit averaged light 
curves in the 15-50 keV energy band from this coded mask in- 
strument are available on a dedicated website 



The Gamma-ray Burst Monitor (GBM; |von Kienlin et aL 



2004 Meegan et al.||2007J > onboard Fermi observes the sky in 
the hard X-ray and soft y-ray regime (about 8 keV to ~30 MeV). 
It permanently provides complete coverage of the unocculted 
sky. Due to its strongly limited spatial resolution, the bright- 
ness of individual sources cannot be determined directly and the 
Earth occultation method is applied ( Case et aL]|201 1| [Wilson- 1 



Hodge et al. 20121. In this work we use the publicly available 
quick look Fermi GBM Earth occultation results provided by 
the Fermi GBM Earth occultation Guest Investigation teams at 
NASA/MSFC and LSL0 which consist of light curves with a 
1 d resolution in four energy bands between 12 keV and 300 keV 
(12-25 keV, 25-50 keV, 50-100keV, and 1 00-300 keV) starting 
from MJD 54690. On average each measurement of CygX-1 is 
based on 18 occultations. 

As these data are pre-screened by the respective instru- 
ment teams, no further selection criteria were applied and we 
use all data available from the start of the each mission un- 
til MJD 56240, resulting in 36454 individual measurements for 
BAT, 9794 for MAXI, and 1443 for GBM (Fig.[l}. 

3. Identifying the states of Cygnus X-1 

3. 1 . General source behavior 

The different periods of source activity and therefore the differ- 
ent population of the individual source states strongly affect our 
ability to distinguish the regions the respective states occupy on 
hardness-intensity diagrams (HIDs) or in other spaces. As a first 
step for our analysis we therefore present an overview of the 
light curves used for this work and over the ASM 5-12 keV to 
1.5-3 keV hardness (Fig. |T]>. 

The source behavior from early 1996 (RXTE launch) until 
the end of 2004 has been discussed by Wil ms et al.| (|2006 ), who 
used a crude definition of ASM based states using the count 
rate only. This classification is sufficient to distinguish main 
activity patterns, and also consistent with more detailed studies 
( |Zdzia rski et a l.|2011| l. We extend this earlier work and supple- 
ment the ASM data with BAT, MAXI, and GBM. By eye, i.e., 
without using the state definitions introduced later in this section 
and shown in color on Fig.[T] we are able to distinguish five main 
periods with different source activity patterns: 

- a pronounced soft state episode in 1996 (up to ~MJD 50350, 
period i), 

- two and a half years of a mainly stable hard state with only 
short softenings, seen as spiky features on Fig.[T[ between the 
end of 1996 and early 1998 (~MJD 50350-51000, period n), 

- a series of failed state transitions and soft states which started 
in early 1998 (Pottschmi dt et al.| 2003) and included a pro- 
longed, very soft period between the end 2001 and the end 
of 2002. This activity period continued until mid-2006 
(-MJD 51000-53900, period in), 



3 http : //maxi . riken . jp/top/index . php?cid=l&jname= 
J1958+352 



** http : //swift . gsfc . nasa . gov/docs/swift/results/ 
|transients/| 



http : //heastro . phys . lsu . edu/ gbm/ 
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- an almost continuous hard state from mid-2006 to mid-2010 
which includes the hardest spectral states ever observed in 
CygX-1 (see |Nowaketll^[20TT1 ~MJD 53900-55375, pe- 
riod iv), 

- a series of prolonged soft states which followed an abrupt 
state transition at ~MJD 55375 and continued during the 
writing of this paper (>MJD 55375, period v). 

Periods i to m are fully and period iv is almost fully covered by 
the ASM. Some ASM data exist during period v, but are affected 
by problems described in Sect. 2.1 and therefore excluded from 



our analysis. 

During period iv the ASM hardness shows higher variability 
than in previous hard states, but the same effect is also seen in 
the hardness of the Crab and is instrumental (Pig. |2j. 

BAT coverage started at the end of period m on MJD 53414 
(mid-February 2005) and continues during the writing of this 
paper. Thus, BAT data are available during soft states, namely 
in period v, but simultaneous coverage of soft states by BAT and 
ASM is lacking. 

GBM and MAXI coverage started MJD 54690 (2008 Au- 
gust 12) and MJD 55058 (2009 August 15), respectively, and 
was ongoing as of time of writing. These instruments mainly 
cover the final phase of period iv and the ongoing period v. No 
simultaneous coverage with ASM exists during the intermediate 
and soft states, except during the phase of ASM deterioration. 

During the joint coverage by ASM and BAT, CygX-1 dis- 
played two softening episodes, which did not reach a stable 
soft state (~MJD 53800-53900 and around MJD 55000). Both 
episodes are clearly visible in the ASM band but not in the BAT 
light curve. This is a remarkable contrast to the full state tran- 
sition of ~MJD 55375, which is associated with a clear drop in 
the BAT count rate simultaneous with an increase of ASM count 
rate. Similar softenings visible in the soft X-ray band but not 
in the hard band have been observed in 1997-1999 by Zdziarski 
et al.| (|2002[), who discuss changes seen in the ASM where no 



corresponding changes were observed in the 20-300 keV band 
with BATSE. Only further long term observations in both soft 
and hard X-rays can help to decide whether only successful state 
transitions are accompanied by a change in hard X-ray flux and 
whether the behavior of the hard component above 15-20keV 
can be used to decide whether a state transition will fail or not. 

After MJD 55200 reliable ASM data are lacking. The BAT 
light curve suggests that while the source left the soft state for 
short periods of time, those hard periods were softer than the 
prolonged hard states in periods n and iv. 

A striking feature of the hardness during the hard states are 
values which exceed the average by a factor of ~5 and more 
( jBalucihska-Church et aLl|2"0"00"l [Poutanen et al.||2008) . They 
correspond to so-called X-ray dips, where blobs of cold material 
in the line of sight cover the source. Since it is virtually impos- 
sible to decide whether an individual ASM dwell was measured 
during a dip, we do not treat these data separately for the purpose 
of this paper (but see |Boroson & Vrtilek|2010| l. 

3.2. r, defined states 



Wil ms et al.| ( |2006[ ) have shown that the spectral shape of 
Cyg X- 1 as observed by PC A and HEXTE (from here on: RXTE 
spectrum) can be well described by an empirical model which 
consists of a broken power law modified by a high energy cut- 
off, absorption described by the tbnew model, see, e.g., Hanke 
[etaLl ( |2"00"9"l i, and the abundances of |Wilms etaLl ( |2000| >. The 
iron Ka line is described by an additive Gaussian emission line 



Table 1. Spearman's rank correlation coefficients p for the correlations 
between the ASM counts and the soft photon index Ti of fits to pointed 
RXTE observations for the cases T[ < 2.7 and T t > 2.7 in ASM energy 
bands 



Energy band [keV] 


1.5-12 


1.5-3.0 


3.0-5.0 


5.0-12 


P (Ti < 2.7) 


0.77 


0.81 


0.76 


0.62 


P (Ti > 2.7) 


-0.60 


-0.54 


-0.55 


-0.68 



at ~6.4keV. The intermediate and soft states additionally show 
a soft excess, which is usually described by accretion disk emis- 
sion ( |Mitsuda et al.|1984[[Makishima et al.|1986| >. The soft pho- 
ton index, Ti, of the broken power law shows strong correlations 
with other spectral and timing parameters on timescales from 
hours to years across the whole range of its values and is a good 
proxy for the source state (e.g., Pottschmid t et al.|2003 Wilms 



et al. 2006; Bo ck et al.|201 1[ |Nowak et al.|201 1) 

We model all spectra both with and without a disk and accept 
the disk as real if the addition of the disk component improves 
the^ 2 value by more than 5"%^] With this approach, we are able 
to achieve good fits (^ 2 ed < 1 .2) for almost all spectra, with a few 
outliers not exceeding x^ ed — 2.5. Less than 18% of the spectra 
with Y x < 2.0 require a disk, but 97% of the fits with F] > 2.0 
do. This agrees with the known behavior of the disk in black 
hole binaries in the different states (e.g., Beiloni|2010 l. For all 
best fit models, we obtain Fi > F2 and £break ~ 10keV. For a 
more detailed discussion of the fits to then available data and 



examples of typical spectra we refer to Wilms et al. (2006). 

Figure [3] shows the dependence of the total rms between 
0.125 and 256 Hz in both high and low energy bands and the 
time lag (averaged between 3.2 and 10 Hz) on the soft photon 
index F\. The timing and spectral behavior is clearly related, but 
complex with changes in both slope and sign. For clarity of the 
plot we do not show error bars; the uncertainity of single mea- 
surements is, however, on the order of or smaller (in the case 
of the hard observations: much smaller) than the spread of the 
correlation at any given frequency. Fj ~ 2, where spectral mod- 
els with a disk become dominant, corresponds to a bend in the 
timing correlations. Another clear kink can be seen at Ti ~ 2.5. 

We use both our spectral and timing analysis results to de- 
fine Ti ranges for different source states that are characterized 
by spectral and timing behavior which is similar within a state 
but different between the three states. Hard states correspond to 
Fi < 2.0, intermediate states correspond to 2.0 < F\ < 2.5, and 
soft states correspond to Fi > 2.5. 

3.3. Simultaneous ASM mapping 

For direct classification of the ASM data we use the 1424 in- 
stances when ASM observations of Cyg X-l are simultaneous 
with spectra from pointed RXTE observations, i.e., fall within 
the good time intervals used for the orbit-wise spectral ex- 
tractions. This corresponds to 2400 individual ASM measure- 
ments since the source is often observed by more than one SSC. 
We treat the measurements of different SSCs independently for 
the actual mapping since instrument alignment onboard RXTE 
would otherwise introduce a bias towards a higher number of 
ASM measurements with Cyg X-l in the field of view of two 
SSCs during pointed observations of Cyg X-l. 



6 Since the PCA data in the lower channels are dominated by system- 
atic errors for a source as bright as Cyg X-l, it is not possible to adopt 
a significance-based criterion for the improvement in^ 2 . 
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Fig. 3. X-ray time lag and fractional RMS as a function of the soft 
photon index T t . 



Figure [4] shows the dependence of the ASM count rate from 
individual ASM-SSC dwells on the soft photon index Ti of the 
broken power law fits. Note that due to the use of individual 
ASM SSC dwells, an RXTE spectrum and therefore a Tj -value 
can correspond to more than one simultaneous ASM measure- 
ment. While the ASM count rate in all bands is correlated with 
Ti, this correlation is not unique: it appears to change sign at 
Ti ~ 2.7 (Table [TJ. Consistent with this observation, Zdziarski 
|et al.| ( |20TT| ) note that on several instances in the soft state their 
derived bolometric flux is lower than during the high-flux hard 
stat^J Thus, cuts in ASM count rate alone cannot separate the 
states. 



Zdziarski et al. 



(2011 



define the soft state according to a power- 
law photon index deduced from the ASM count rates. This approach 
does not exactly correspond to the approach chosen here, however, it is 
similar enough to allow rough comparisons. 
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Fig. 4. Dependence of the ASM count rate in different bands on 
the soft photon index of simultaneous pointed RXTE observations. 
Shown are total ASM counts (1.5-12keV, black dots) and counts in the 
three ASM energy bands: band A (1.5-3 keV, orange upside down trian- 
gles), band B (3-5 keV, green triangles) and band C (5-12 keV, magenta 
squares). 



The Ti-based definition of states introduced in Sect. |3.2| re- 
sults in in 1608 ASM measurements in the H -defined hard state, 
in 455 ASM measurements in the Fi -defined intermediate state, 
and in 337 ASM measurements in the Fi -defined soft state. Each 
of the states primarily populates a well-defined distinct area in 
the ASM HID (Fig. [5]). A visual inspection reveals that nei- 
ther cuts in hardness only nor cuts in count rate only yield a 
good division between the states: the softest observations tend 
to have count rates in the transitional state range (~40-50cps), 
the hard state area extends to low hardnesses usually associated 
with transitional or even soft states, and the hard and transitional 
state strongly overlap both in ASM hardness and count rate. To 
account for these features we choose the following ansatz: an 
ASM observation with count rate c and (5-12 ke V/l. 5-3 keV) 
hardness h is defined as hard if c ^ 20 cps regardless of the hard- 
ness. For c > 20 cps the cuts between the states are defined by 
linear functions of the form 



casm = '«asm,/ • (h - ho,d 



(1) 



where i e {hard, soft). The line with slope mASM,hard an d x ~ 
intersection Zio.hard divides the hard and the intermediate state, 
and the line with slope w so ft and intersection fto, S oft divides the 
intermediate and the soft state, respectively. 

We determine the best division between the states such that 
the fractional contamination of the ASM-defined states by dif- 
ferent Ti-defined states is minimized. Contamination here is de- 
fined for the hard state (and accordingly for soft and transitional 
states) as the fraction of all measurements classified as hard us- 
ing the ASM that is classified as transitional or soft according 
to T\ . Initial fits indicate that good separations of the states are 
achieved for fto,hard ~ ^o,soft- We therefore reduce the number of 
free parameters for the cuts and set 



h 



O.hard 



h 



O.soft 



(2) 
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Table 2. Overview over all sky monitor based state definitions for Cyg X- 1 . 



State 


ASM-based a 




MAXI-based b BAT-based c 


GBM-based d 


hard 

interm. 

soft 


c < 20 V c < 55 ■ 
c> 20 A 55 ■ (h - 
c> 20 A c> 350 


(h - h ) 

-h )<c< 350 ■ (h - h ) 
■ (h - h ) 


c M <l-4-/i M 

1.4-fiM < cm < 8/3 -hu ■■■ e 

8/3 ■ /im < cm cb < 0.09 counts crrr 2 s _1 


f 
f 

/ < 0.6 Crab 



Notes. <a) With ASM 1.5-12keV count rate c in counts s -1 , ASM (5-12 ke V/l. 5-3 keV) hardness h, and/i = 0.28. (b) With MAXI2^keV count 
rate c M in counts s _I and MAXI (4-10 keV/2-4 keV) hardness h M . (c) With BAT normalized 15-50 keV countrate c B in counts cirr 2 s _I . <d) With 
daily GBM 25-50 keV flux /. (c) Discriminating between hard and intermediate states is not possible based on BAT lightcurves alone. The source 
is defined as in BAT-based hard or intermediate state for c B > 0.09 counts crrr 2 s~'. (f) Discriminating between hard and intermediate states is not 
possible based on GBM daily lightcurves alone. The source is defined as in GBM-based hard or intermediate state for / > 0.6 Crab. 
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Fig. 5. PCA to ASM mapping. Grey data points are all ASM mea- 
surements of Cyg X-l in the shown range. Blue points represent PCA 
defined hard states, green intermediate and red soft states. Black lines 
show the cuts defining the states in the ASM HID. The light blue shaded 
region corresponds to the position of the hard state in the HID, the light 
red shaded to the position of the soft state. The intermediate state region 
is shown without shading. 
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Fig. 6. Percentage of contamination of the ASM defined states (hard 
state shown in blue, intermediate in green, and soft in red) for time 
delays At = ±(0-1.5) h, ±(1.5-3) h, ±(3-4.5) h, etc., between the ASM 
and the pointed RXTE measurements. Simple linear fits to the data are 
shown as dashed lines to illustrate the overall trends. 



For the best cuts we obtain ho = 0.28, mh ar d = 55, and m SO ft = 
350 (Fig.|5]andTable|2). 

The contamination by other redefined states is <5% for the 
hard state, <10% for the intermediate state and <3% for the soft 
state. Note that since the source behavior changes continuously 
from one state to the other, such that the states do not represent 
three distinct, fully independent regimes, we do not expect a per- 
fect separation of the states (e.g., Wilms et al. 2006 ). The spread 
in count rate and hardness is amplified by the orbital and super- 



orbital modulations of Cyg X-l and by dips (e.g., Poutanen et al. 
2008). The stronger contamination of the intermediate state is 
expected: as a transitional state between the hard and the soft 
state, it is short-lived and confined by two divisional lines. The 
separation between the hard and the intermediate state appears 
especially unclean; we therefore advise to treat the classification 
cautiously when an observation is close to this cut. 

To test our approach we compare the ASM-based behaviour 
with the results of the spectro-timing analysis of the quick, 
observationally exceptionally well covered intermediate to soft 
transition presented by Boc k et al.| ( |2011| >. In particular, we 



can recover the moment of the transitions at slightly before 
MJD 53410. 



3.4. Non-simultaneous ASM mapping 

In general, an all sky monitor and a pointed instrument will not 
observe a source at the same time, so that we need to assess how 
well a given monitor pointing can be used to characterize the 
source state during a non-simultaneous observation. 

For every RXTE spectrum we consider ASM measurements 
within 1.5 h intervals At = +(0-1 .5) h, ±(1.5-3) h, ±(3-4.5) h, 
etc., up to 48 hours. The length of the intervals is motivated 
by the length of the RXTE orbit of ~1.5 h, during about half of 
which Cyg X-l is visible. For simultaneous ASM measurements 
with different SSCs we use the average for all following analysis. 
We obtain 135901 pairs of ASM and RXTE state classifications 
(the same ASM measurement may be used to classify several 
RXTE spectra, if the RXTE spectra are close enough). For every 
delay interval and for every ASM defined state, we determine the 
percentage of spectra with a different RXTE classification. 
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Fig. 7. PCA to MAXI mapping. Gray data points are all MAXI 
orbitwise measurements of CygX-1 in the shown energy range. Blue 
circles represent PCA defined hard states, green intermediate states and 
red soft states. Black lines define the state cuts. The light blue shaded 
region corresponds to the hard state in the MAXI HID, the light red 
shaded region to the soft state and the region without shading to the 
intermediate state. 



Figure[6]shows that this contamination remains stable for the 
hard state for up to ~48 h. For the soft state, the contamination 
reaches 20% for a 48 h delay and is even larger for the interme- 
diate state, which also shows a very strong increase within the 
first ~10h. Note that strictly simultaneous data, which were dis- 



cussed in Sect. 3.3 are not taken into account here, resulting in 
larger starting contaminations. The results are similar when us- 
ing positive or negative delays only. The trends are expected, 
because the hard state occurred often in long, stable stretches 
during the RXTE lifetime, while the intermediate state is short 
lived due to its transitional nature. 



3.5. MAXI mapping 

Since only 37 MAXI measurements are strictly simultaneous 
with pointed RXTE data, we use MAXI data within At = 1 h be- 
fore and after a pointed RXTE observation and obtain 219 MAXI 
measurements with H -based state classifications, which offer a 
better overall statistics with 107 hard states, 101 soft states, and 
1 1 intermediate states. Increasing At to 2 h does not yield a bet- 
ter coverage of the intermediate state, neither does attempting 
to map ASM onto MAXI, since a simultaneous coverage by the 
two instruments is only available during the end of period iv. 

To find the best approach to define MAXI based states we 
consider all possible different combinations of count rates mea- 
sured with MAXI (the three energy bands as introduced in 
Sect. |2.3| and the overall count rate) and different combinations 
of hardness measures within the MAXI bands. The clearest sep- 
aration can be achieved when using the ratio between count rates 
in the medium (4-10 keV) and the low (2-4 keV) MAXI bands 
(MAXI-hardness, /imaxi) and the count rate in the low (2—4 keV) 
MAXI band, c M axi, as done in Fig. [7] The MAXI 4-10keV/2- 
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Fig. 8. Histogram of BAT fluxes of T| -defined states: the hard state is 
shown in blue, the intermediate state in green, and the soft state in red. 



4 keV hardness is also the closest correspondence to ASM hard- 
ness we can achieve using publicly available MAXI light curves. 

The sparse coverage of the intermediate state makes it hard 
to separate the three basic states. Given that the hard and the soft 
state populate distinct parts of the MAXI HID (see Fig. [7]l and 
knowing the shape of the cuts in the ASM HID, we separate the 
states by two linear functions of the form 



CMAXI - 7WMAXI,; " ^MAXI, 



(3) 



where i e {hard, soft) and where mhard separates the hard and 
the intermediate state and m so f t separates the intermediate and 
the soft state. The absence of an x-intersection and of a thresh- 
old count rate value is motivated by the lower number of classi- 
fied data points compared to the ASM HID, which also prevents 
us from direct fits for m. The best values obtained by eye are 
«MAXi,hard = 1-4 and m M Axi,soft = 8/3 (Table |2j. 

For a reader interested in performing her own classification 
using MAXI data, we note that these are conservative cuts in the 
sense that we obtain the purest intermediate state possible here. 
As in the case of RXTE-ASM we expect the separation between 
the hard and intermediate state to be especially unclean. 

3.6. BAT mapping 

Swift-BAT light curves in different energy bands are not read- 
ily available. Since BAT also does not cover soft X-ray energies 
where the contribution of the accretion disk becomes important 
as the spectrum softens, we use only the light curves in the stan- 
dard 15-50keV band for our analysis and therefore do not per- 
form a two-dimensional mapping. 

Since only 290 BAT measurements are strictly simultane- 
ous with pointed RXTE data, we employ the same approach as 
for MAXI (Sect. [33} and use BAT data within At = 1 h before 
and after the RXTE spectrum to obtain 1339 BAT measurements 
with an RXTE based state classification: 895 in the hard state, 
124 in the intermediate state, and 320 in the soft state. 

Only 290 BAT measurements are strictly simultaneous with 
pointed RXTE data. Following our study of the reliability 
of non-simultaneous data (Sect. |3.4| i, we use BAT data within 
At = 1 h before and after a pointed RXTE observation and ob- 
tain 1339 BAT measurements with T\ -based state classifications, 
which offer a better overall statistic. 

Figure [8] shows the histogram of the BAT fluxes in the dif- 
ferent states. While the soft and hard state show clearly different 
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BAT fluxes, the intermediate states populate the same region as 
the hard ones. The BAT measurements, therefore, do not enable 
us to separate the hard and the intermediate states. The soft state 
can still be identified by a cut at a BAT area normalized count 
rate of cbat = 0.09 counts cm -2 s _1 For BAT count rates below 
cbat CygX-1 is in the soft state, for fluxes above it in the hard 
or intermediate state (Table [2]). The contamination by the other 
state is ~7% in both cases. 

We also considered an ASM to BAT mapping where we clas- 
sified BAT data by using the closest ASM measurements if one 
exists within ±0.5 h around the BAT measurement in order to 
increase the number of data points with classification. No simul- 
taneous good ASM and BAT data are available in the soft state 
(see Fig. [JJand Sect. 2.1 1. The histograms of the hard and in- 
termediate state follow the trends apparent from Fig. [Sj it is not 
possible to separate the two states using BAT flux. 



3.7. GBM mapping 

Since GBM lightcurves are publicly available in daily bins only, 
strictly simultaneous mapping between RXTE spectra and GBM 
is not possible. We therefore define an RXTE spectrum as simul- 
taneous to a GBM measurement if the good time interval used 
to extract the spectrum lies within the GBM measurement. If all 
simultaneous RXTE spectra show the same state, the GBM mea- 
surement is classified as belonging to this state. If the states of 
the RXTE spectra differ, the GBM measurement is not classified. 

The 104 GBM measurements with simultaneous RXTE ob- 
servations include on average 3.8 (reaching from 1 to 8) indi- 
vidual RXTE spectra. Only 6 out of the 104 measurements are 
unclassified due to state ambiguities. The remaining valid clas- 
sifications include 61 hard states, 3 intermediate states, and 34 
soft states. All 6 unclassified measurements include hard as well 
as intermediate state classifications. There were no GBM mea- 
surements including both intermediate and soft, or hard and soft 
RXTE spectra. The distribution of these results is consistent with 
CygX-1 spending most of the time covered by GBM observa- 
tions in stable hard and soft states. That 6 out of 9 GBM mea- 
surements including a -defined intermediate state also include 
Ti -defined hard states reflects the high variability and instability 
of the intermediate states. 

We first consider GBM HIDs. These show a clear separa- 
tion in two regions: a region with mainly soft and a region with 
hard and intermediate observations, which are clearly separated 
in GBM flux only. We therefore consider GBM fluxes only and 
calculate histograms of GBM fluxes of Ti -defined states. The 
best separation between the states can be achieved when using 
25-50 keV fluxes (Fig. [9]). While hard and intermediate states 
cannot be separated, the soft state can be divided from the other 
two by a cut at a GBM flux of 0.6 Crab: for fluxes <0.6 Crab the 
source is in the soft state, for fluxes >0.6 Crab in hard or inter- 
mediate state (Table |2j. We do not give values for contamination 
here because of the long integration times of the GBM and be- 
cause we excluded states where we know the source to change 
the Ti defined state during one GBM measurement. We also 
note, that since the hard and intermediate states cannot be dis- 
tinguished using GBM, we do not introduce additional bias into 
our classification by not including GBM-measurements with si- 
multaneous RXTE spectra in both hard and intermediate states. 



4. The states of Cyg X-1 

4.1. The statistics of CygX-1 states 

Using the classification from the different instruments derived in 
the previous sections, we can assess the activity pattern during 
the different periods defined in Sect. |3.1| A detailed break-down 
of the occurence of different states is given in Table [3] 

Period i contained a prolonged soft state and the source 
spent twice as much time in soft states than in the intermediate 
state during this time. We lack all sky monitor coverage before 
MJD 50087 and therefore do not know whether this soft period 
observed was a part of a longer series of soft states. 

In period m the source spent as much time in the intermediate 
state as in the soft state, in accordance with the observation of 
multiple failed state transitions (e.g., |Pottschmidt et al.|2003] ). 

Periods n and iv showed similar activity patterns dominated 
by a long, stable hard state, although data from the ASM were 
already affected by deterioration of the instrument in period iv 
(Sect. [XT). The state classifications based on the ASM and the 
BAT agree well for period iv. 

The disagreement between the statistics derived from Swift- 
BAT and MAXI classification for period v first seems worri- 
some. A visual inspection of Fig. [TJ indicates, however, that 
this mismatch may be due to the gaps in the MAXI light curve. 
Indeed, out of the 7794 BAT measurements which fall into 
MAXI gaps (defined here as times when the interval between 
two MAXI measurements is larger than 6 h), 3876 are hard or 
intermediate and 3918 are soft, while 60% of all BAT measure- 
ments during period v are soft. MAXI gaps therefore fall more 
often onto hard/intermediate states than onto soft states during 
this period and the difference seen in Table [3] is due to the in- 
complete MAXI coverage. MAXI data can therefore be used to 
determine states at certain times, but not to investigate the over- 
all frequency of the occurrence of different states. A similar cau- 
tion also applies to the GBM. Here, the disagreement between 
the GBM data and other instruments is partly due to the use of 
daily lightcurves, where the short term variability is averaged 
out (Fig. [TJ. This can be directly seen from the higher resolution 
BAT data. Using individual BAT points during period v, 60% 
of all measurements are classified as coming from the soft state 
(Table BJ. Classifying the daily BAT fluxes, however, increases 
the fraction of soft states to 68%, which is consistent with the 
GBM considering the uncertainties in both the GBM- and BAT- 
based definitions (Figs. [8] and [9). Data with a time resolution 
<1 d are therefore crucial for a reliable state classification. 

Summarizing these observations, we estimate the typical 
inter-instrument systematic error of the state determination to be 
better than 10% for the post-RXTE instruments. 



4.2. Stability of States 



Sect. |4. 1 1 highlights that the different states are clearly different 
in their stability. While the source tends to remain in the hard 
and the soft state for a prolonged time, not unexpectedly the in- 
termediate state is far more unstable. This well known general 
behavior can be quantified using the RXTE-ASM data, the data 
set in our measurements that covers the longest time span, does 
not have gaps, and can be reliably used to define all three states. 

To calculate the state statistics, for each ASM data point 
measured at a reference time, f ref , we calculate the probability 
that Cyg X- 1 has changed its state after a certain interval, and the 
probability that the source remained in the same state throughout 
this interval. More formally, we calculate the transited fraction, 
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Fig. 9. Left: Daily GBM HID. S denotes the 12-25 keV flux, H the 25-50 keV flux. The gray dashed line denotes the zero flux level. Total 
GBM measurements are shown in gray. Right: Histogram of daily GBM fluxes of Ti -defined states. The hard state is shown in blue, intermediate 
state in green, and soft state in red. Total GBM measurement are shown in black. The gray dashed line denotes the zero flux level. 

Table 3. Time CygX-1 spent in the different states (hard - H, intermediate - 1, soft - S) as measured by the RXTE-ASM, Swift-BAT, MAXI and 
Fermi-GBM all sky monitors. 



Period 2 


MJD 


RXTE-ASM 




MAXI 




Swift-BAT 


Fermi-GBM 






H 


I 


S 


H 


I 


S 


H&I 


S 


H&I S 


I 


50087 b -50350 


37% 


19% 


44% 














II 


50350-51000 


99% 


1.0% 


0% 














III 


51000-53900 


63% 


20% 


17% 














IV 


53900-55375 


97% d 


2% d 


<l% d 








99% 


1% 




V 


55375-56240 c 








17% e 


8% e 


75 % e 


40% 


60% 


23% f 77% f 



Notes. Between MJD 50087 and 55200 CygX-1 spent a total of 75.6% of its time in the hard state, 12.8% in the intermediate state and 11.6% in 
the soft state as measured with ASM. 

Percentages are only shown if at least most of the period is covered by the respective instrument. {a) Periods defined in Sect. 3.1 (b) Start of period i 
defined by the start of the RXTE-ASM measurements. (c) The end of period v is defined by the availability of data at the time of writing. <d) ASM 
data only until MJD 55200. (c) MAXI data affected by gaps in the lightcurve. (f) Daily average values are used for GBM. 



^trans, i-e., the probability that an observation made a certain time 
after the reference observation will find CygX-1 in a different 



source state (Fig. 10 1. The probability that the source state is un- 
changed is 1 - /Vans- We determine P t i- a ns by calculating the frac- 
tion of all measurements in a different state than the reference 
measurement made at times f; with Afi < tj - f re f < Af 2 where 
{Afi,Af 2 } e {{Oh, 1.5h},{1.5h,3.0h},{3.0h,4.5h},...}. This ap- 
proach is equivalent to that employed in Sect. |3.4| where we used 
RXTE spectra as reference measurements. 

Ptrans does not take into account the possibility that the source 
might have undergone state changes between f ref and f,. To ad- 
dress this possibility we calculate the cumulative transited frac- 
tion, P ct ,„, which measures the probability that at least one source 
transition has happened up to a time interval At after the refer- 
ence measurement. In order to assess the influence of possible 
outliers and/or misclassifications, we define a state transition by 
the existence of more than n e {0, . . . , 12} measurements in a dif- 
ferent state than the reference measurement in the time interval 
f re f < t < f re f + At. The probability that the source has remained 
in the same state for at least a time Af is then given by 1 - P ct „ 

(Fig. [nj. 

As expected from a simple look at the lightcurves (see Fig.[TJi 
hard states are the most stable states: in more than 90% of 
all cases, 200 h after a given hard state ASM measurement the 
source is still in the hard state (Fig. [TO}. In over 30% of all cases 
a period of 300 h after any given hard state measurement will 



not contain even a single measurement in a different state. This 
number increases to over 70% for larger n, i.e., when very short 
excursions to the soft/intermediate state and/or outliers due to 



misclassifications are ignored (Fig. 1 1 



Though less stable and shorter than hard states, soft states 
generally show a similar behavior. Even for n = 0, a significant 
number of soft state measurements is not followed by a state 
transition within 300 h, i.e., prolonged stable, soft states lasting 
more than 12 days exist. Such states are observed by ASM in pe- 
riods i and m (Fig. [TJ. MAXI and BAT data suggest an increased 
occurrence of such stable soft states also in period v. 

The intermediate state behaves differently. The transited 
fraction strongly increases within the first few hours. Akeady 
fifty hours after a reference intermediate measurement about half 
of the measurement will show the source in a different state. The 
cumulative transited fraction approaches 100% even for n — 12 
at 300 d (Fig. [TT[ >. Intermediate states are therefore short and 
unstable compared to hard and soft states. 



5. Summary 

Based on pointed RXTE observations, we provided criteria to 
define X-ray states (hard, intermediate, and soft) using light 
curves from all sky monitor instruments. In particular we have 
shown that: 
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At [h] 

Fig. 10. Probability, Ptransi that the source state is found changed 
a time interval At after a previous state determination for the different 
states (blue: hard, green: intermediate, red: soft) using ASM data. 




At [h] 

Fig. 11. The probability, 1 - P ct „, that the source has remained in 
the same state (blue: hard, green: intermediate, red: soft) for at least 
the time At for different values of n of ignored number of possible out- 
liers/misclassifications (n = 0: solid, n = 4: dashed, n = 8: dot-dashed, 
n = 12: dot-dot-dashed). The gray dashed line represents 0. 

- due to the complex source behavior, simple state definitions 
based just on the source count rate or just on the hardness do 
not adequately describe states, 

- a combination of RXTE-ASM total count rate and hardness 
can be used to define states before MJD 55200 (see Table [2] 
for exact cut), 

- the best separation of states is achieved with simultaneous 
ASM data. Data within At ±6 h result in a contamination of 



<10% for the hard and soft state and <30% for the interme- 
diate state, 

- a combination of MAXI count rate and hardness can be used 
to define states (see Table[2]for exact cuts), 

- soft coverage is necessary to define states and the lack of 
such does not allow to distinguish the hard and the interme- 
diate state using publicly available BAT light curves. BAT 
light curves can, however, be used to distinguish the soft 
state from hard/intermediate states: the source is in the soft 
state below a threshold of 0.09 counts cirT 2 s _1 , and in the 
hard/intermediate state above it (Table [2}, 

- the lack of soft coverage does not allow a separation be- 
tween hard and intermediate state with publicly available 
GBM data, with the analysis being further hindered by the 
low time resolution of the GBM lightcurve, but a rough 
state classification is to consider CygX-1 in the soft state 
for fluxes <0.6Crab and in the hard/intermediate state above 
that threshold (Table [2}, 

- the hard state is by far the most stable state of Cyg X-l, 
followed by the soft state. The probability that the source 
remains in the hard state for at least one week (200 h) is 
>85% (using P c t,i2)- Soft states are slightly less stable, but 
the probability of a soft state being longer than one week is 
still ~75%. Intermediate states are short-lived and typically 
last a few days at most, implying that they can only be caught 
with monitor data with a time resolution of better than 1 d. 

The state classification introduced here can be reliably used to 
define source states where no X-ray continuum spectrum mea- 
surement is available. The high frequency of all sky monitor 
measurements for RXTE-ASM, Swift-BAT, and MAXI enables 
us to catch short flares and especially quick state transitions. 
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